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The deuteron breakup on heavy target has been investigated in the frame of an improved quantum molec- 
ular dynamics model, focusing on the production of neutrons near zero degree. The experimental differential 
cross sections of neutron production in 102 MeV d+C reactions are reproduced by simulations. Based on the 
consistency between the model prediction and experiment, the feasibility of producing neutron beam through 
the breakup of deuteron on carbon target has been demonstrated. Because of the nucleon Fermi motion inside 
the deuteron, the energy spectrum of the inclusive neutron near 0° in laboratory exhibits a considerable energy 
broadening in the main peak, while the long tail at low energy side is suppressed. By measuring coincidentally 
the accompanying proton from deuteron breakup, the energy of the neutron can be tagged with an intrinsic un- 
certainty of about 5% (1c). It enables the application of well-defined energy neutron beam in an event-by-event 


scheme. 


PACS numbers: 29.25.Dz, 25.45.-z, 24.10.-i 


Introduction - High-quality monochromatic neutron beam 
has very important applications in many fields. For instance, 
the neutron radiation hardness evaluation for the electronic el- 
ements in aviation devices requires neutron beam with well 
defined energy [1, 2]. The assessment of the neutron dose 
received by the patients in heavy ion therapy or by the astro- 
nauts in the aircraft relies on the calibration using well de- 
fined neutron beam [3]. In fundamental researches, high ener- 
gy neutron beam is of significance for the calibration of high 
energy neutron detector, as well as for accumulating data of 
the cross section of neutron-induced reaction [4]. Undoubt- 
edly, high-quality and high-energy neutron beam becomes a 
necessary tool [5-8]. Since neutron has no charge, while the 
neutron beam below 20 MeV can be produced by nuclear re- 
action with fixed Q value at given angle, monochromatic neu- 
tron beam above 20 MeV is usually produced in laboratory 
with proton induced reaction on Lithium target [7, 9, 10]. S- 
ince the nucleons are bounded in the target nuclei via nuclear 
force, the knocked-out neutron possesses a long energy tail 
which can not be eliminated [7, 11], nor be tagged based on 
event-by-event scheme. 


The properties of deuteron nucleus and its applications have 
been the subject of intense studies since its discovery in 1932 
[12]. One year later, deuteron beam was produced for the 
first time [13]. The first photodisintegrations of deuteron was 
measured [14]. Using photographic photometry, the spin of 
deuteron was determined as J = 1 [15]. Deuteron was found 
to possess an electric quadrupole moment [16], indicating the 
existence of additional D wave component. Such D wave 
could be generated by the tensor part of the one-pion exchange 
(OPE) potential [17] and led to profound understanding of the 
nature of nuclear force. Since 1950s, enormous experiments 
of ed scattering have been performed to extract the form fac- 
tors of deuteron nucleus, including charge monopole, magnet- 
ic dipole and charge quadrupole form factors. For a review, 


one refers to [18]. 

Deuteron-induced reactions exhibits some intriguing fea- 
tures too. Coulomb polarization of deuteron-induced trans- 
mutation was discovered in 1935 [19, 20]. Such polarization 
is reported to increase the fusion reaction rate very recently 
[21]. The measurement of deuteron breakup can be dated back 
in 1970s [22]. Moreover, The isovector reorientation effec- 
t has been predicted in deuteron-induced scattering on heavy 
target, providing a novel probe to nuclear symmetry energy of 
hot discussions so far [23, 24]. 

Can deuteron be used to generate high-quality neutron 
beam? The answer seems to be yes. It was mentioned that 
deuteron induced reaction can be used to produce neutron 
beam and the neutron energy can be tagged in the specific 
p(d, np)p channel [18], and some early measurements on the 
n — p angular and energy correlation in deuteron breakup has 
been reported [25]. Recently, Jin et al. reported that deuteron- 
induced spallation reactions could be used as a powerful way 
to generate neutron beam [26]. These results are due to the 
fact that deuteron is a loosely bound nucleus with binding 
energy Eg = 2.2 MeV and is easily disassociated in the re- 
action with target. Our motivation in this paper is to study 
the breakup process of deuteron in peripheral reaction, focus- 
ing on the advantageous possibility to precisely determine the 
energy of the neutron beam. Since monochromatic deuteron 
beam is easily to gain on modern accelerator, if the idea to 
produce neutron beam with deuteron is feasible, it is expected 
that the neutron keeps the half of momentum of the incident 
deuteron and can be used as a (quasi)monochromatic neutron 
beam. The paper is organized as following. Section 2 presents 
briefly the improved quantum molecular dynamics (ImQMD) 
model. Section 3 discusses the distribution properties of the 
generated neutron beam. Section 4 is the conclusion. 

Model Description - InQMD model is an extend version 
of QMD[27]. A detailed description of the ImQMD model 


and its version ImQMD05 and their applications can be found 
in Refs. [23, 28-32], where deuteron- and nucleon- induced 
reactions are particularly simulated. Within the ImQMD05 
model, the nucleon is represented by a Gaussian wave packet 
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where r; and p; are the centers of wave packet of the i" nu- 
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where p, Pn, Pp are the nucleon, neutron, and proton densi- 
ty, 6 = (Pn — Pp)/(Pn + Pp) is the isospin asymmetry degree. 
The parameters in Eq. (2), except C, and y which describe 
symmetry potential energy, are fully determined by Skyrme 
interactions. The Skyrme parameter set MSLO [33], one of 
Skyrme parameter sets which satisfy the current understand- 
ing of the physics of nuclear matter over a wide range of ap- 
plications [34], is used in the calculations. The parameters are 
listed in Table I. 


While the initialization of the heavy target nuclei is done as 
usual as that in traditional QMD, the deuteron is initialized se- 
mi classically in a simplified scheme as that in [23]. At the end 
of the ImQMD calculations, clusters are recognized by a min- 
imum spanning tree (MST) algorithm [27, 35] widely used in 
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where Dyin = 0.5 fm, bmax = 8 fm and Ab = 0.5 fm are tak- 
en in our calculation. The possibility distribution is written as 
SE, Q, bi) = y (E, Q, bi) /N where y is the yield of the neutron 
in each cell and N is the total simulated events under a certain 
reaction condition. It is shown on the scattering plot that there 
are two components standing out in statistics. One distributes 
nearly homogeneously in wide angular range with E, < 20 
MeV, originating from the target region, the other situates in 
projectile region keeping approximately the beam energy per 
nucleon at very forward angle with ab < 10°. Except for 
the two components, the neutrons originating from the target 
fragmentation are distributed widely on the plot. The intensity 
in the target rapidity region is not high, this is because i) we 
chose a light target containing less neutrons and ii) the bind- 
ing energy of deuteron is much less than the separation energy 
of neutron in carbon nucleus. 


cleon in the coordinate and momentum space, respectively. 
According to the results from our previous studies, a tradition 
value of 02 = 2.0 fm? is appropriate for the intermediate en- 
ergy nucleon-induced reactions. Nucleons in a system move 
under the mean-field with the nuclear potential energy density 
functional, which reads 


8/3 C, p! 


7 + E (Woy + EE Von =P) + br ag + 6", (2) 
Po Po Po 


the QMD calculations. In this work, the nucleons with rela- 
tive momenta smaller than 250 MeV/c and relative distances 
smaller than 3.0 fm are coalesced into the same cluster. The 
information of cluster with excited energy is input into the s- 
tatistical decay model GEMINI[36, 37] to perform statistical 
decay calculations. 

Results and Discussions - We first calculate the phase space 
distribution of free neutrons in d+C reactions at 102 MeV 
where the experimental data are available[38]. We choose the 
stable carbon instead of heavy metal as the target simply for 
the purpose to suppress the contributions of the neutrons from 
target. Fig. 1 presents the doubly differential cross section of 
do-/dQdE, as a function of the kinetic energy En and polar an- 
gle @ in laboratory for the neutrons produced in the reactions. 
The differential cross section is calculated by 
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The energy spectra at a certain angle in laboratory can be 
obtained from Fig. 1 by slicing the 2-dimensional histogram. 
Figure 2 presents the neutron energy spectra with logarithmic 
scale on the ordinate at various angles in comparison to the 
experimental data [38]. It is shown that our calculation re- 
produces well the main peak in the vicinity of beam energy 
per nucleon, particularly at small angles with ia» < 5°. The 
height of the main peak decreases rapidly with 6), in accor- 
dance with the data, indicating that the neutrons are peaked at 
forward angle. On the other hand, due to the clustering inef- 
ficiency of transport model, i.e., the model counts less (more) 
clusters (nucleons) than experiment, the yield of neutron at 
very low energy, originating mainly from the target fragmen- 
tation, are overestimated by approximately 30%. Besides, a 
very small peak near 100 MeV, which is the direct knock-out 
neutron from the target, is not reproduced in our calculation. 


TABLE I. Parameter set used in the ImQMD calculations. 
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FIG. 1. (Color online) The doubly differential cross section of 
do/dQdE, as a function of the kinetic energy and emitted angle with 
beam direction 


If one count only the main peak of the spectrum at 0° and 5°, a 
full width at half maximum (FWHM) of approximate 20 MeV 
can be derived. It indicates readily that the direct breakup of 
deuteron can be used as a source of neutron beam with the 
broadening of 20 MeV in energy. 

Figure 3 further presents the neutron energy spectra in the 
angular range 6, < 2° produced in d+C reactions with the 
incident deuteron energies of 102, 200 and 400 MeV, respec- 
tively. The results of the ImQMD simulation with a GEMINI 
afterburner is represented by the histograms in the upper pan- 
el. The standard deviation og, of the neutron beam normal- 
ized to the average neutron energy (En) is plotted in panel (b). 
At Ebeam = 200 MeV, the oz, reflecting the monochromaticity 
of the secondly neutron beam is about 33 MeV, which is larg- 
er than the neutron spectra produced in p+Li reactions[7, 11]. 
The long tail at the low energy side is much suppressed in d+C 
compared to that in p+Li channel of neutron production. 

Although the binding energy of deuteron is only 2.2 MeV, 
the proton and neutron in the deuteron have Fermi motion with 
an average momentum of 75 MeV/c, which causes the broad- 
ening of the neutron energy in the main peak. In the breakup 
or stripping reaction of deuteron, the total energy of the pro- 
jectile keeps approximately constant, anti-correlation is then 
expected between the energy of the proton and the neutron 
from deuteron breakup. Figure 4 (a) presents the scattering 
plot of the proton energy vs. the neutron energy in d+C reac- 
tions at various beam energies ranging from 60 to 400 MeV. It 
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FIG. 2. (Color online) Comparison of neutron spectra at various an- 
gles between experimental results [38] and ImQMD calculations. 


is evident that a sharp band is presented at each energy point, 
confirming the anti-correlation of the proton and neutron en- 
ergy. The total energy of the proton and the neutron, En + Ep, 
is displayed correspondingly in panel (b), where the neutrons 
are counted within lab < 2°. It is evident the total energy is 
constant with an intrinsic FWHM about 5 MeV despite of a 
tiny tail appears at low energy side. This feature has been ob- 
served in the d+!?C at 270 MeV where the correlated protons 
and neutrons are measured at 0° [40]. According to the exper- 
iment, the FWHM of the En + Ep spectrum at 0° is about 2.5 
MeV. Moreover, the intrinsic broadening on the total energy 
shows insignificant dependence on the beam energy, implying 
that it originates from the Fermi motion of the nucleons rather 
than from the kinetic effect. It suggests that if one can mea- 
sure the accompanying proton at forward angle, the energy 
of the neutron can be tagged with greatly improved accuracy, 
compared to the raw main peak in Fig. 3. 

A further question remains that whether one can main- 
tain high efficiency of delivering a neutron beam with well- 
determined energy using the proton tag. Fig. 5 (a) and 
(b) presents the efficiency of the proton tag Rag and the 
monochromaticity of the neutron beam using the coincident 
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FIG. 3. (Color online) The energy spectra (a) of neutron in d+C at 
Eq = 102, 200 and 400 MeV deuteron beam energy within Oa, < 
2°. The relative standard deviation o,,/E, as a function of neutron 
energy is plotted in panel (b). 


proton in different angule cut as a tag, respectively. Here Rag 
is derived by the ratio of the number of coincidence between 
p and n over the total number of neutrons. When a neutron is 
found in flab < 2°, one searches the coincident proton with an 
angular cut in laboratory as indicated in the figure. It reads 


Poxp 


Riag = P 
n 


(4) 

It is shown that the efficiency increases rapidly at low 
beam energy and gradually saturates at 90% with beam en- 
ergy above 200 MeV/u. Comparing the results at different 6, 
cut implies that the coincident proton is still at forward angles 
in accordance with the experimental observation[25]. The 
monochromaticity in panel (b) is defined as the standard devi- 
ation C£, of the total energy En + Ep normalized to the mean 
neutron energy, TE! (Ea). Clearly to see, the monochro- 
maticity increases with the beam energy and is better than 3% 
in the production of neutron beam above 200 MeV. 

Finally, Fig. 6 compares the cross section of the neutron 
production in two different channels, the deuteron-induced 
production d+!°C in this work and the conventional channel 
of p+'Li [11, 42-44]. The angular emittance of the neutron 
beam are both 6, < 2° with respect to the primary deuteron 
beam in laboratory system. It is shown that, the cross section 
of neutron production in d+!?C is higher by a factor of 100, 
slightly dependent on the energy. This result has importan- 
t implication, i.e., in order to obtain the same neutron beam, 
the primary beam intensity of deuteron with (d,np) channel 
can be lower than that of proton with (p,n) channel by 2 order 
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FIG. 4. (Color online) The correlation of high energetic neutron and 
proton in the forward angle (a), the total energy distribution of the 
neutron and the accompanying proton (b). 


of magnitude. As an advantage, the background is reduced 
significantly. The date point of En = 51 MeV is sitting on 
top of the simulation[38], indicating our calculation is reli- 
able. It suggests that the breakup of deuteron provides a novel 
method to generate high energy neutron beam of well deter- 
mined energy and efficiency. It is worth mentioning that the 
cross section of the secondary neutron from deuteron breakup 
increases with the atomic number of the target [40], but likely 
in the price of producing more neutrons originated from the 
target, which degrade the monochromaticity. It shall be noted 
here that in our transport model calculations, the mechanisms 
of stripping and breakup are not distinguished, since no cou- 
pled channels are considered [41]. 


Summary - To summarize, we calculated the neutron spec- 
tra produced in d+C reactions at various incident energies in 
comparison with experimental data if available. It is shown 
that the neutron originating from the breakup of deuteron 
keeps the kinetic memory of the projectile with some broad- 
ening on the energy arising from the Fermi motion of the nu- 
cleon in the projectile. It is shown that the accompanying pro- 
ton from the deuteron breakup exhibits sharp anti-correlation 
with the neutron and the total energy of the neutron and the 
proton keeps constant with less than 5 MeV intrinsic energy 
variation. Thus, with an efficiency of about 90%, the energy 
of the neutron can be tagged by the accompanying proton on 
an event-by-event basis with such enhanced precision. The 
cross section of producing the neutron beam is much higher 
in (d,np) channel than in (p,n) channel. 
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FIG. 5. (Color online) The efficiency of proton tag and the 
monochromaticity of total energy of neutron and proton. The neu- 
trons in 6, are counted, while the coincident protons are searched in 
different angular cut on 6, in laboratory. 
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FIG. 6. (Color online) The cross section of the neutron production 
with two channels d +!” C and p +’ Li. 
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